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ABSTRACT 



We use Galaxy Zoo 2 visual classifications to study the morphological signatures 
of interaction between similar-mass galaxy pairs in the Sloan Digital Sky Survey. 
We find that many observable features correlate with projected pair separation; not 
only obvious indicators of merging, disturbance and tidal tails, but also more regular 
features, such as spiral arms and bars. These trends are robustly quantified, using a 
control sample to account for observational biases, producing measurements of the 
strength and separation scale of various morphological responses to pair interaction. 
For example, we find that the presence of spiral features is enhanced at scales < 70 
hjQ kpc, probably due to both increased star formation and the formation of tidal 
tails. On the other hand, the likelihood of identifying a bar decreases significantly in 
pairs with separations < 30 hj^ kpc, suggesting that bars are suppressed by close 
interactions between galaxies of similar mass. 

We go on to show how morphological indicators of physical interactions provide 
a way of significantly refining standard estimates for the frequency of close pair in- 
teractions, based on velocity offset and projected separation. The presence of loosely 
wound spiral arms is found to be a particularly reliable signal of an interaction, for 
projected pair separations up to ~ 100 hj^ kpc. We use this indicator to demonstrate 
our method, constraining the fraction of low-redshift galaxies in truly interacting pairs, 
with M* > 1O 9 ' 5 M0 and mass ratio < 4, to be between 0.4 - 2.7 per cent. 

Key words: galaxies: general — galaxies: evolution — galaxies: interactions — 
galaxies: structure 
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1 INTRODUCTION 

When galaxies approach one another, their mutual grav- 
itational attraction can result in substantial disruptions 
to their morphologies, such as tidal arms, counter arms, 
bridges and tails. Many studies have shown, both an- 
alytically and using numerical simulations, that galax- 
ies of similar mass can provoke dramatic disturbances 
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in the stellar distributions of one another, with the de- 
tails depend ing on the orbital parameters of the inter- 
actio n (e.g., iToomre & Toomrd 19721; IBarnes &: Hernquist 



19921 : iHoward et all Il993l : iGerber fe Lambl Il994l : lBar~ 
20 111 ). Gravitational perturbations can also redistribute 



the gas content of galaxies, potentially lea ding to c h anges 
in their star- format ion properties (e.g., iNoguchil Il988l ; 
IBarnes fe Hernauistl [l996). These effects are strong func- 
tions of pair separation, and hence should be most obvi- 
ous after the galaxies' first pass, and particularly around 
times of closest approach. Many interactions are ultimately 
likely to result in the complete morphological transformation 
of th e galaxies involved (e.g., lToomrelll977l ; iHopkins et al.l 
120081 ). However, more subtle effects are expected both ear- 
lier, while the galaxies are on their initial appr oach, and at 
times of wide separat i on between passes (e.g., Perez et all 

201 ih . 



times or wide s eparation betwee n passes (e.g., 
120061 : lLotz et al.|[2005 IStruck. Dobbs. fc Hwang 



The expected strength and prevalence of pair interac- 
tions mean they are potentially important for determining 
the properties and evolution of the galaxy population. It is 
therefore critical that we test our theoretical expectations 
of the effects of such interactions by studying representative 
samples of interacting systems. Furthermore, we may utilise 
the observed frequency of interacting pairs to discriminate 
between the details of cosmological galaxy formation mod- 
els. 

A relatively straightforward, and physically motivated, 
definition of 'interacting' galaxies is a pair for which the tidal 
force experienced across one of the galaxies, averaged over 
its internal dynamical time, Ft, is at least some specified 
fraction of the gravitational force binding the outer regions 
of that galaxy, F g . By this definition, all interacting pairs 
should produce significant internal dynamical effects, which 
would have otherwise been absent, in at least one of the 
member galaxies. In practice, however, it is difficult to mea- 
sure the forces involved. We may estimate them by studying 
the effects of an interaction, but the relative orientations and 
types of the galaxies in each pair, as well as observational 
limitations, lead to large variations in the apparent effects 
for interactions of a given strength, Ft/Fg. 

A more convenient definition, which is roughly equiva- 
lent, although only statistically applicable, is that a pair are 
'interacting' if their gravitational influence upon one another 
could have observable effects in a favourable orientation and 
mix of galaxy types. For example, a pair of elliptical galax- 
ies might not display signatures of an interaction in a given 
observation, but would still count as 'interacting' by this def- 
inition if the tidal forces they are experiencing would have 
been sufficient to produce an observable signature in a pair 
of spirals. The observational details of a particular dataset 
therefore fix the minimum Ft/F g probed. This definition re- 
moves much of the incompleteness associated with only con- 
sidering pairs with observational signs of interaction, but of 
course only a fraction of such 'interacting' galaxies will pos- 
sess observational signatures. As we shall see later in this 
paper, it is nevertheless possible to constrain the fraction of 
galaxies interacting according to this definition. 

Studies of galaxy pairs typically discuss close pairs, 
bound pairs, merging pairs, or pairs with observational dis- 
turbances, but often mix the usage and definitions of these 
classes. Interacting pairs are closely related to bound pairs, 
for which the sum of the gravitational potential energy and 



kinetic energy of the pair is negative. However, not all in- 
teracting pairs are bound, particularly where both are part 
of a larger system, such as a galaxy cluster. Likewise, not 
all bound pairs will be experiencing significant tidal inter- 
actions. Interacting pairs are also closely related to mergers. 
Galaxies in bound pairs will typically merge on relatively 
short time-scales if they are experiencing significant tidal 
interactions, as the kinetic energy of the pair orbit is trans- 



ferred to deiormmg and heatin g the i: 
tion of each galaxy (e.g. Struck 1999) 



We can identify close pairs of galaxies which are likely 
to be sufficiently near to one another such that they 
are interacting (and potentially bound and will eventually 
merge) using projected di stance and line-of-sight velocity 
([Charlton &; Salpeterlll99lh . However, this approach suffers 
from significant contamination and incompleteness (with re- 
spect to the above definitions of physically meaningful in- 
teracting, bound or merging selections), due to a lack of full 
spatial information and the inverse relationship between rel- 
ative velocity and separation for loose pairs (i.e. very close 
pairs can have very large relative peculiar velocities, so ap- 
pear significantly separated in redshift space) . Observational 
signatures of interactions, for example visual classifications, 
quantitative morphological measurements or induced star 
formation, may be used to improve the selection of truly in- 
teracting galaxies. However, as often one wishes to study the 
physical effect of interactions, one must be careful to avoid 
a circular argument. 



The early atlases o f lVorontsov-Velvaminovl (1 19591 . 1 19771 ) 
and lArpj (1966) clearly demonstrated that interactions be- 
tween galaxies can have profound effects on their morpholo- 
gies, providing examples of bridges, tails, distorted spiral 
patterns, and other features. These morphological changes 
where observed between pairs up to over 100 kpc, as 
is the case for the bridging fila ment of Arp 295. The re - 
stricted three-body simulations of IToomre &; Toomrd (|l972l ) 
clearly demonstrated that these features are the result of 
strong tidal forces between the interacting galaxies. That 
galaxy intera ctions can also induce st ar- format ion was first 
suggested bv lLarson fe Tinslevl (|l978l ). who found that the 
scatter in the UBV colours of interacting galaxies from the 
Arp atlas was sign ificantly larger than normal galaxies from 
the Hubble Atlas (|Sandagd[l96ll ). Similar evidence for in- 
teraction induced star- format ion has also been found over 
a wide r ange of the energy spectrum, from near- UV to ra- 
dio ( e g. IKennicutt fc Keell 19841: iKeel et alJll985b iBushouse 



1990) 



19861: IKennicutt et al.lll98 7: Bushouse 198 <1: iHummel et al 



Studies which use quantitative measures of morphol- 
ogy find signs of interaction at fairly small proj e cted sep- 
arations. Using the CAS method of IConselicel (|2003l ) to 
measure galaxy asymmetry (A) and concentration (C), 
iHernandez- Toledo et al.l (|2005l ) found that both these quan- 
tities increase, relative to isolated galaxies, for galaxy pairs 
with sep arations less than the p hotometric diameter of the 
primary. |Pe Propris et al.l ((2007) reliably identified interact- 
ing pairs with projected separations up to r p < 50 kpc 
using A > 0.35 and visual confirmation, for a sample of 
pairs with line-of-sight velocity differences AV < 500 km 
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s _1 Q Similarly, lEllison etal show that asymmetry 

increases for r p < 50 kpc for a sample of pairs with 
AV < 200 km s" 1 . 

Meanwhile, studies which probe the effects of tidal inter- 
action through star formation modulations find ch anges up 
to larger projected separations. iNikolic et al.l (|2004 ) demon- 
strate an increase in star formation at at r p < 70 Ji^q kpc for 
early and mixed type pairs, and at r p < 430 /if x kpc (their 
maximum separation probed) for late type pairs. They also 
find that pairs with r p < 110 kpc show a strong increase 
in central con centration, suggestive of nuclear starbursts. 
iLi et al.l (|2QQ8h find a star formation increase for close pairs, 
with star formation rate (SFR) enhanced by a factor of 1.5 
at r p < 140 h^o kpc to a factor of 4 at r p < 30 kpc. 
This strong dependence on r p , is contrasted with a weak 
dependence on luminosity ratio, with the star formation 
enhancement being strongest in lower lu mino sity galaxies. 
lEllison et all (|2008h jRobaina et al.1 (|2009h and lPatton et al.1 
(|2011) all find a strong in crease in SFR for r p < 40 h^Q kpc, 
while Patton et al ] (l201lh also see a smaller increase up to at 



least r p < 80 /if kpc (their maximum separation probed) 
for AV < 200 km s _1 pairs. There is also evidence that 
equivalent levels of tidally induced st ar formation require 
smaller r p in dense r environments (e.g. jLambas et al.ll2003l : 
lAlonso et alll2004l ). 

Some of the variation in the separation scale at which 
different studies begin to see the effects of tidal interactions 
is likely due to differences in the mass and luminosity ranges 
of the samples, as well as the methods used. However, it ap- 
pears clear that the effects of tidal interactions are found at 
larger projected separations when identified by induced star 
formation (up to 80 to 100 h^o kpc), compared to quantita- 
tive measurements of asymmetry (up to ~ 30 kp c) . This 
is consistent with the results of lLotz et al.l (| 20081 ). which 
used simulations to show that quantitative morphological 
methods for finding merging galaxies, such as A, Gini and 
M20, are most sensitive for galaxies undergoing close pas- 
sages and during the post merger phase. Induced star for- 
mation, on the other hand, will be evident between passes, 
when the galaxies achieve a wide separation before falling 
back towards one another, or in galaxies which have experi- 
enced a close encounter but will not merge. Note that in dry 
mergers there may be no star- format ion signature of the in- 
teraction, and morphologic al features will t ypically only be 
observable for short times ([Bell et al.ll2006h . 

As mentioned previously, interacting galaxies can pro- 
duce distinctive morphological features such as tidal arms, 
counter arms, bridges and tails, which are best classified vi- 
sually. These features are extremely reliable for discriminat- 
ing truly interacting galaxies from interlopers in close pair 
catalogues. Features such as two loosely- wound tidal arms 
may not be detectable by quantitative morphology meth- 
ods because these galaxies may not appear to be sufficiently 
asymmetric or disturbed, especially between the first and 
second pass when the galaxies may appear to be widely 
separated. One of the advantages of using visual morpho- 
logical classifications over automated methods is the abil- 
ity to identify very faint and subtle features. Tidal features 



1 Throughout this introduction, separations quoted from other 
studies have been converted to units of kpc as necessary. 



are known to become rapidly undetectab le as a function 
of time and survey im aging depth fe.g., [Bell et al.l [2006; 
ISchawinski et al.ll2010h . however we find that the Galaxy 
Zoo classifications used in this paper are extremely sensitive 
to faint features. Furthermore, as we will show, by study- 
ing the occurrence of such features in a statistical sense, 
and making weak assumptions concerning the observability 
of physical interactions, we are able to make decisive state- 
ments concerning the prevalence of interactions. 

The visual classification of peculiar, disturbed and in- 
teract ing galaxies has a long history, beginning withlHubblel 
(Il926l ) . The catalogu es of Vorontsov- Velvaminovl J 1 959 ' 



Il97/t ). ISandagd (Il96lh . and lArpl (ll966T ) complied together 
a significant number of galaxies with obvious tidal fea- 
tur es. Such work continues to be valua ble today, for exam- 
ple iBridge, Carlberg. fc Sullivanl (|201(J ) uses visual classifi- 
cations of galactic bridges and tails in the CFHT Legacy 
Survey to study the evolution of the galaxy interaction frac- 
tion (GIF) with redshift. These galaxies were either isolated 
merger remnants or fairly close interacting pairs, due to their 
requirement tha t galax y pa irs be connected by a bridge. 
iNakamura et al.1 (I2003T ) and iFukugita et~ai1 (I2007T ) visually 
classified a subsample of ~ 2500 bright galaxies from the 
Sloan Digital Sky Survey ( SPSS: I York et alJl2000h imaging 
of SDSS galaxy objects, finding that ~ 1.5 per cent of galax- 
ies in their nearby magnitude-limi ted sample show morpho - 
logical indications of interaction. iNair & Abraham! (|2010h 
provide an impressive catalogue of detailed visual classifica- 
tions for over 14 000 bright SDSS galaxies, which includes 
information regarding tidal tails and other indicators of in- 
teraction. 

The Galaxy Zoo project has enabled visual classifica- 
tion to be performed for extremely large samples, allowing 
the cont inued use of t his va luable technique with modern 
surveys. ISkibba et all (|2009h obtained the marked correla- 
tion function for the Galaxy Zoo 1 merger classification like- 
lihood and found that it increases sharply in their closest r p 
bin (of 170 h^Q kpc width), and found evidence that most 
of this increase was for pairs with r p < 30 kpc. Taking 
an alter native approach. iDarg et al.l (|2010ah and lDarg et al.l 
(2010b) imposed thresholds to select Galaxy Zoo 1 galaxies 
with high merger classification likelihoods and study their 
frequency and properties. Most of these galaxies are either 
highly disturbed systems or very close pairs. While these 
studies have been successful at identifying a subset of in- 
teracting pairs, they primarily select galaxies which have 
relatively small projected separations and so do not typi- 
cally identify interacting pairs which are at large projected 
separation between their first and second close passes. 

In this paper we use visual classifications from Galaxy 
Zoo 2 to study what morphological changes are taking place 
in interacting galaxy pairs as a function of physical projected 
separation (r p ) and line of sight velocity difference (AV). 
We then use these results to estimate the frequency of pair 
galaxy interactions in the local universe. 

In Section [2] we describe the data set and sample se- 
lection, in Section [3] we outline the methods employed, in 
Section [4] we present our results, and in Section [5] we sum- 
marise our results and discuss their implications. A AC DM 
cosmology is assumed throughout, with Qa = 0.7, Q m = 0.3, 
and h 70 = H o /(70 km s" 1 Mpc" 1 ). 
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2 GALAXY ZOO 2 DATA AND SAMPLE 
SELECTION 

Although there have been many attempts at completely au- 
tomating morphological classification, visual inspection re- 
mains the preferred method for many astronomers. How- 
ever, for the large samples produced by modern surveys, vi- 
sual classification is not feasible for a normal research team 
to pe rform in a reasonable time. Galaxy Zoo ([Lintott et al.l 
l20Q8h is an online citizen science project designed to ad- 
dress this problem, by involving large numbers of the pub- 
lic in classifying the morphological features of galaxies. 
The original Galaxy Zoo website collected classifications for 
nearly one million galaxies from SPSS Data Release 6 (DR6; 
lAdelman-McCarthv et al.l l2008). which were processed to 
produce catalogues of visual classifications together with es- 
timates of their accuracy. These data have been relea sed to 
the publicQ, and are described in lLintott et al.l (|201lh . 

The original Galaxy Zoo (GZ1) was limited to coarse 
morphological classification. Following its success, a subse- 
quent project was launched, this time collecting much more 
detailed morphological information via a question tree, for 
a subset of ~ 300 000 of the brightest galaxies from Galaxy 
Zoo. This project, named Galaxy Zoo 2 (GZ2), ran from 16 
February 2009 until 22 April 2010 and collected 16 340 298 
classifications (comprising a total of 58 719 719 questions) by 
83 943 participants for 325 651 galaxy imageqj. Both GZ1 
and GZ2 used gri composite colour images pr ovided by the 
SPSS , created following the prescription of iLupton et al.l 
(|2004h . These were disp layed such that each galaxy had the 
same apparent size. See lMasters et all ([201 lh for additional 
discussion of the GZ2 dataset. 

Figure [T] presents the classification tree used for GZ2, 
including the actual button images shown to the GZ2 par- 
ticipants. These images attempt to symbolise the answer 
to each question. Most of the participants have no formal 
astronomy training and, although there is a tutorial which 
presents real examples for each answer, it is likely that they 
rely on these images to a reasonable degree while making 
their classifications. 

The GZ2 sample selection includes 273 783 galaxies 
from SDSS DR7 with Rqo > 3" and rp e tro,AB < 17.0 and 
located in the contiguous North Galactic Cap regiorQ- In 
this paper we require redshift information, and therefore 
restrict our sample to galaxies with spectra in DR7. The 
overall spectros copic complete n ess of this GZ2 sample is 86 
per cent. As in iMasters et al.l (|2010h , we limit our sample 
of galaxies to z < 0.09 to ensure sufficient resolution for 
spiral classifications and thereby red uce somewhat the red - 
shift classification bias discussed in Ba mford et al.l (|2009). 
A lower redshift limit of z > 0.01 is imposed to avoid the 
de-blending of large nearby objects into multiple photomet- 
ric objects by the SDSS pipeline, and ensure that redshift- 
derived distance moduli are sufficiently accurate. Applying 
these limits results in a sample of 148 291 galaxies. 



2 Galaxy Zoo data is publicly available at 
http: / /data.galaxyzoo.org/ 

6 The GZ2 website is archived at http://zoo2.galaxyzoo.org 
4 GZ2 also included ~30,000 galaxies from both the normal- and 
coadd-depth imaging of the SDSS Stripe 82 region, but these data 
are not used in this paper. 



We make use of the MPA- JHU DR7 median stellar mass 
measurement^!, which are based on fits to the SPSS pho - 
tom etry, following the methods of iKauffmann et al.l |2003) 
and lSalim et all (|2007T ). 



3 METHOD 

3.1 Pair and control samples 

In this study we focus on galaxy pairs with stellar mass 
ratios between 1:1 and 4:1. These would amount to major 
mergers, if the interaction were to proceed that far. The 
candidate galaxy pairs in our parent sample are selected to 
have an absolute line of sight velocity difference of AV < 
5000 km s _1 and a projected physical separation (converted 
from an angular separation using the average redshift of the 
pair) of r p < 1000 kpc. 

To probe the effects of changing AV between pairs we 
divide our sample into 100 km s _1 bins for AV < 500 km 
s _1 , as well as two larger AV bins with 500 < AV < 1000 
km s" 1 and 1000 < AV < 5000 km s" 1 . The 1000 < AV < 
5000 km s _1 pairs serve as our control sample, as they should 
all be physically unassociated. These are used to account 
for biases that result in a dependence of some morpholog- 
ical classifications on the projected separation of galaxies, 
irrespective of any true interaction. They also indicate the 
typical level of morphological features present in the general 
galaxy population, for comparison with our pair samples. 

Note that pairs with AV > 5000 km s _1 have a sig- 
nificantly different luminosity distribution (skewed to more 
luminous galaxies) compared to the lower AV samples. 
This is a result of AV approaching the velocity range of 
the survey, and hence the sample becoming dominated by 
more distant, intrinsically brighter galaxies. We judge the 
1000 < AV < 5000 km s _1 range to be a good compro- 
mise between maintaining consistent luminosity selections, 
while maximising the control sample size and minimising its 
contamination by physically associated pairs. 

3.2 Morphology probabilities 

GZ2 participants are asked a series of questions for each 
image, with each answer determining the subsequent ques- 
tion, as depicted in Fig. [1] We record the individual answers 
provided and, after some down- weighting of the most incon- 
sistent participants, use these to construct a catalogue. This 
catalogue gives the total number of times each galaxy was 
presented and, for each question, the fraction of votes for 
each possible answer, which we denote f(A). The f(A) for 
all A belonging to a single question must, of course, sum 
to unity. To give a concrete example, consider that a given 
galaxy has been presented to 20 participants, 16 of which an- 
swered 'Features or disk' to the initial question 'Is the galaxy 
simply smooth and rounded, with no sign of a disk?'. This 
galaxy thus has a vote fraction of f (Features) — 0.8. The 16 
participants who answered Features were then asked 'Could 
this be a disk viewed edge-on?', with 12 answering 'No', and 
thus f(Edge-on=No) = 0.75. All these 12 were subsequently 

5 this catalogue is publicly available from 
http:/ /www. mpa-garching.mpg.de/SDSS/DR7/Data/stellarmass. html 
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Figure 1. The Galaxy Zoo 2 question tree 



asked 'Is there a sign of a bar feature through the centre of 
the galaxy?', resulting in a split of / (Bar — Yes) =0.5 and 
f(Bar=No) = 0.5. Again, all 12 were then asked 'Is there 
any sign of a spiral arm pattern?', with a 'Yes' answer frac- 
tion of f(Spiral=Yes) = 0.75. The question tree continues, 
and we have only considered the route through the ques- 
tion tree taken by the majority, but this is sufficient for the 
reader to fully understand both the traversal of the question 
tree and the definition of f(A). 

The vote fraction f(A) may, very roughly, be consid- 
ered to represent probabilities regarding a galaxy's mor- 
phology. The uncertainty expressed by these probabilities 
results from a combination of limited observational informa- 
tion, true morphologies which do not precisely align with the 
possible answers, and differing interpretations by each par- 
ticipant in judging the correspondence between the image 
and each answer. In this case f(A) — p(A\M), the condi- 
tional probability of the galaxy having morphological feature 
A given that it possesses the set of morphological features 
iVf , since each participant must have identified the galaxy 
as having the preceding features in the question tree, M, in 
order to be asked about feature A. 

The conditional probability compares the likelihood of a 
galaxy having feature A against the alternative answers for a 
single question. However, it does not necessarily give a good 
representation of the presence of a particular morphological 
feature. An object may have a high f(A), but still be unlikely 
to possess morphological attribute A. To assess the overall 
likelihood of A, one may calculate an estimate for the joint 



probability of A and M. Formally this is, 



p(A n M) 



■ p(A\M)p(M) 

■ f(A) n £/( a )' 

QCM aeQ 



(i) 

(2) 



where Q are subsets of the set of answers M partitioned by 
question, and the a are the individual answers in Q. 

As we do not distinguish between the different paths 
which lead to a given question (though this would be pos- 
sible to do from the raw data), we only consider cases for 
which M represents the sum of all possible paths to an- 
swer A. In this case, we denote p(A) — p(A D M), where 
p(A) is the probability of morphological feature A together 
with any M for which asking the question with answer A is 
relevant, as denned by the question tree in Fig. [1] Remem- 
ber that these probabilities only include the observational 
information available in GZ2, so do not necessarily equate 
to the true probability of a particular morphology feature. 
For example, for a galaxy which has f (Features = Yes) < 1 
or f(Edge-on=No) < 1, the value of p(Spiral=Yes) does not 
take into account the unobservable (in terms of GZ2) pres- 
ence of spiral arms in apparently smooth or edge-on galaxies. 

For the example above, 



p(Spiral=Yes | Features D Edge-on=No) 
f(Spiral=Yes) 



0.75, 



(3) 
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whereas, 

p(Spiral=Yes) = 
p(Spiral=Yes fl Features D Edge-on=No) = 
f(Spiral=Yes)f(Edge-on=No)f (Features) = 0.45. 



(4) 



This indicates that there is only a moderate probability that 
the object in question has visible spiral arms, although if 
one is willing to accept that it does have features and is not 
edge-on, then it probably does possess spiral arms. Notice 
that the Bar question has been implicitly omitted from this 
calculation as all its possible answers continue on to the 
Spiral question, and hence it would contribute a factor of 
[f(Bar=Yes) + f(Bar=No)] = 1. 

Whether one works with f(A) or p(A) depends upon 
the question one is considering. It is particularly useful to 
examine trends in terms of f(A) itself, as this reflects the 
behaviour of a specific morphological feature, irrespective of 
other morphological variations. However, due to the limited 
total number of times each object is considered, when p(A) 
is low, f(A) will be highly quantised and subject to high 
Poisson noise. In this case, one can consider only objects for 
which f(A) is meaningful, by imposing a minimum threshold 
on p(M), which we denote by f(A\p(M) > £), for some 
threshold t. 

To study the dependence of f(A) on projected sepa- 
ration (r p ) for a particular sample, we take the mean vote 
fraction f(A) in each bin of r p , 



f(r P ,A) = ±jrf(A) 



(5) 



where the sums are over the N galaxies in each r p bin for 
that sample. This is done for each of the answers, A, shown 
in Fig. [U 

For many answers we find that the control sample dis- 
plays a dependence on r p . Given their velocity separation, 
these trends cannot be due to any physical interactions 
within the pairs. Instead they must arise from the appar- 
ent close projection of the galaxies, and may be considered 
to be a 'projection bias', which will contaminate any signa- 
ture of physical interaction in the lower AV pairs. In order 
to remove this contamination, we remove the control sample 
trends versus r p from the observed trends for the other pair 
samples. If we regard the f(r p ,A) as conditional probabili- 
ties, the control sample trend may be removed by, 



(/p-/c)/(l-/c): 



(6) 



where F is the conditional probability in the absence of 
the projection bias, f p is the measured conditional prob- 
ability for a sample of physically associated pairs, and f c 
is the measured conditional probability for the control sam- 
ple. F(r p ,A) therefore represents the conditional probabil- 
ity, P(r p , A\M), of morphological feature A being observed 
as a result of the galaxy being in a pair with separation r p , 
given that the galaxy displays morphological features M. 
The uncertainties on F(r p , A) are given by, 



2 (1 ~ fp) 2 , 

°F = — a fc + 



1 



2 

\2 a fp 



/c)4 -7c ■ (l_ /c )2 ~/P • ( 7 ) 

where cr/ c and <7/ p are determined from the standard error 
of f c and f p in each r p bin. 

Note that for physically associated pairs, with low AV, 



the projected galaxy density increases as r p decreases, such 
that the number of galaxies in a given r p bin, and hence 
the signal-to-noise of f(r p , A) remains reasonably constant. 
However, for unassociated pairs the projected galaxy den- 
sity is constant as a function of r p , and thus at small r p 
the signal-to- noise of /(r p , A) drops substantially. Unfortu- 
nately, this limits the accuracy of the projection bias correc- 
tion and translates into higher uncertainties in F(r p ,A) at 
low r p . 

In order to quantify the strength and scale of the trends, 
we fit the F(r p , A) with the function, 



F(r p ) = aexp(-r p /&) + c, 



(8) 



where a is the size of the change in F from large to small 
r, b specifies the r scale of the trend, and c accounts for 
possible constant systematic offsets between the physically 
associated pairs and control sample. We find that this simple 
empirical function is able to well represent most of the F(r p ) 
trends we observe. Ideally c should be zero, and in any case, 
for F(r p ) to represent a probability, < c < (1 — a). How- 
ever, systematic offsets of c in either direction are possible 
due to small differences in the sample selections of physically 
selected pairs and the control sample. These are difficult to 
avoid for different AV selections, but should not have an 
r p dependence, justifying the use of a simple constant to ac- 
count for them. Reassuringly, we find that c is generally very 
close to zero, signifying that sample selection differences are 
indeed minimal. The fitting method provides uncertainties 
on a, b and c, which enables us to judge the significance 
of differences in the trends between samples and plot confi- 
dence intervals on the fitting functions. 



3.3 Counting companions 

In Section 14.61 we study the number of close compan- 
ions per galaxy as a function of projected separation, 
jVc(r p ), for a mass- l imited sample. Following the methods of 
IPatton et al . (2000, 2002), weights are applied to account for 
pairs near the survey boundaries (wb 2 ), pairs near the red- 
shift boundaries (w V2 ), global spectroscopic incompleteness 
weights (w s )i as well as angular spectroscopic incomplete- 
ness weights to correct for fibre collisions at small angular 
separation (wg 2 ) . The uncertainties in these weights are de- 
termined primarily from uncertainties in the astrometric and 
redshift measurements. Because the uncertainties in the pair 
statistics obtained in Sect ion H~6l are significantly larger than 
the astrometric and redshift uncertainties, uncertainties are 
not explicitly calculated for these weights. 

In Section T4. 6 1 of this paper we will derive an additional 
weight (itfint) to account for the occurrence of interlopers: 
galaxies in close pairs (as judged by r p and AV), but which 
are not truly interacting. Previous studies typically ignore 
this, or estimate a constant value, whereas we will derive its 
dependence on r p . 

The total weight assigned to each galaxy is thus 



WNo 



W s We 2 Wb 2 W V2 Mint 



(9) 



The total number of companion galaxies of host galaxy 
i, with projected separation r p , is given by summing these 
weights for all companions within a given projected separa- 
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Figure 2. Mean vote fractions for the Odd=Yes answer and subsequent 'odd' categories which show trends versus r p : Merger, Irregular 
and Other. The main plot in each panel shows the mean vote fraction of AV < 500 km s _1 close pairs, corrected for projection bias, i.e. 
F(r p , A) = P(r p , A\M (points with error bars), together with best fit (blue line) to these points using equation [8] and the corresponding 
la confidence region (green lines). The insets plot f(r p , A) for the AV < 500 km s _1 close pairs (blue) and 1000 < AV < 5000 km s _1 
control sample pairs (red), from which F(r p ,A) is calculated. 



tion, 

N 

N c , t (r p ) = ^2R(r p j)w N2 j , (10) 

3 = 1 

where N is the total number of galaxies in the sample, and 
R(r) = 1 if r is in the current r p bin, and otherwise. 

The average number of close companions per galaxy, as 
a function of projected separation, is then calculated as the 
mean of the number of companions of each galaxy, weighted 
by the spectroscopic incompleteness, 

^(r P ) = ^f>c,^p). (11) 

For N c (r p ) <C 1 this is equivalent to the fraction of galaxies 
with close companions. 



4 RESULTS 

Following the method described in Section 13.21 we have 
examined all of the GZ2 answers to ascertain which are 
most relevant to studying galaxy interactions. The ques- 
tions which are of most interest for this paper are those 



regarding odd features, bars, the spiral arm winding tight- 
ness and the number of spiral arms. As can be seen in Fig. [I] 
for each object every participant is asked Ts there anything 
odd?' and, if they answer 'Yes', asked to specify one odd 
feature from seven alternatives. In order for a participant 
to be asked the questions regarding spiral arm number and 
winding tightness they must answer Features, Edge-on=No 
and Spiral— Yes in the preceding set of questions. 

In sections 14.11 14.21 14.31 and 14.41 we use the bias- 
corrected vote fractions, F(r p ,A), corresponding to the 
probability of galaxy property A occurring as a result of 
the galaxy being in a pair with separation r p . For exam- 
ple, F(r p , 1 Arm) = P(r p , 1 Arm \ Features fl Edge- on= No D 
Spiral=Yes). We use the full sample, without applying any 
thresholds for the preceding questions, in order to see what 
classification trends exist for that morphological feature in 
relation to pair separation irrespective of other morpholog- 
ical features. 

In sections 14.51 and 14.61 we select individual galaxies 
based on their morphological features, and wish to mini- 
mize the impact of noisy f(A) values for galaxies with low 
p(A). We therefore impose a threshold on p(A) to select only 
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Figure 3. Mean vote fractions for the GZ1 merger answer. The 
main plot shows the mean vote fraction of AV < 500 km s _1 close 
pairs, corrected for projection bias, i.e. F(r p ,A) = P(r p , A\M 
(points with error bars), together with best fit (blue line) to these 
points using equation [8] and the corresponding la confidence re- 
gion (green lines). The insets plot f(r p , A) for the AV < 500 km 
s _1 close pairs (blue) and 1000 < AV < 5000 km s _1 control 
sample pairs (red), from which F(r p ,A) is calculated. 

objects for which asking the question with answer A is likely 
to be appropriate. 

4.1 The Odd class 

In Fig. [2] we study the trends in 'Odd' GZ2 morpholo- 
gies as a function of pair separation. The main panels plot 
the mean vote fraction, F(r p ,A), of AV < 500 km s _1 
close pairs, correcting for projection bias effects using the 
1000 < AV < 5000 km s _1 control sample, as outlined 
m Section[321 This may be interpreted as the probability, 
P(r p , A\M, of a galaxy in this sample displaying the spec- 
ified morphological attribute, purely as a result of being in 
a pair with separation r p . The raw vote fractions, f(r p , A), 
from which the corrected quantities are determined, are also 
shown in the inset panels. 

Considering first the Odd=Yes answer itself, top left in 
Fig. [21 we see from the inset plot that the raw vote frac- 
tions, f(r p , Odd=Yes) increase strongly with decreasing r p 
for both the AV < 500 km s _1 close pairs and the control 
sample. As the pairs in the control sample are physically 
unassociated, this indicates that GZ2 participants identi- 
fied some galaxies as 'Odd' simply as a result of their ap- 
parent separation. This projection bias strongly contami- 
nates the f(r p , Odd=Yes) results of the AV < 500 km s _1 
close pair sample. However, there are offsets between the 
f(r p , Odd=Yes), such that the Ay < 500 km s _1 close pairs 
are more likely to be marked 'Odd' than the control sample, 
particularly at 20 < r p < 40 Ji^q kpc. 

Indeed, when we correct for the projection bias, we see 
clear evidence for a trend. GZ2 participants were more likely 
to identify an object as 'Odd' if it is in a AV < 500 km s _1 
close pair with small projected separation. The probability 
of a galaxy being labelled 'Odd', as a result of being in a 
pair, P(Odd=Yes), varies from zero for r p > 80 /if x kpc to 
~ 0.5 for the smallest projected separations. The empirical 



function defined in Eqn. [8] does a good job of representing 
this trend. 

If a participant answered Odd= Yes, then they were then 
asked to identify the odd feature more precisely by select- 
ing one option from a variety of possibilities. The remaining 
panels of Fig.[2]plot F(r p , A) = P(r p , A | Odd=Yes) for three 
of these options. From the top right panel of Fig. [2] we see 
that the Merger answer mimics the behaviour of Odd=Yes, 
although note that we are plotting the conditional probabil- 
ity given that the object does display an 'Odd' feature, and 
hence the Odd=Yes behaviour itself is not included in this 
quantity. 

We see, from the raw vote fractions, that the Merger 
answer is strongly dependent on the apparent separation 
of galaxy pairs. Galaxies are often marked as a 'Merger' 
because they appear close together in the image, even when 
there are no other signs of interaction. Despite this, there 
is a clear enhancement of Merger features at 20 < r p < 
40 hjo kpc for AV^ < 500 km s _1 close pairs over that 
seen for the control sample. In the projection bias corrected 
F(r p , Merger) we also see that physically interacting low AV 
pairs do have an additional probability of being identified as 
a merger, although it is a noisy function of r p . 

The Merger answer displays cross talk with the other 
'Odd' categories. As f(r p , Merger) increases with decreas- 
ing r p , the vote fractions of most of the alternative answers 
(i.e. Ring, Arc, Disturbed, Irregular, Dust Lane) decrease, for 
both the low AV and control sample pairs. However, with 
the projection bias accounted for, at the smallest separations 
(r p < 10 h^Q kpc) F(r p , Irregular) and F(r p , Disturbed) in- 
crease, possibly at the expense of F(r p , Merger). This is con- 
sistent with users preferentially classifying separated pairs as 
'Mergers' and interpreting overlapping pairs as single 'Irreg- 
ular' or 'Disturbed' objects. An enhancement of GZ2 Dust 
Lane feat ures in merging early -types has already been dis- 
cussed by IShabala et al.l (|2012h , although we only see a ten- 
tative indication of this with our method. 

Interestingly, the Other category displays a contrasting 
behaviour, showing a decreasing F(r p , Other) with decreas- 
ing r p . This appears to be a result of other, more specific, 
categories being favoured for pairs with a real physical as- 
sociation. Looking at the Other button image in Fig. [1] it is 
clear why some projected close pairs are given this classifi- 
cation, especially non-interacting early type galaxies. 

All these results have important implications for the 
use of the GZ2 'Odd' classifications. In particular, they im- 
ply that for physically interacting close pairs, the likelihood 
of being classified as a merger actually begins to decreases 
at the smallest separations as galaxies begin to be classi- 
fied more as irregular or disturbed. Given the large amount 
of projection bias and the cross-talk between categories, it 
is difficult to use any of the 'Odd' categories alone as in- 
dicators of pair interactions. It seems that using the more 
general Odd=Yes provides the most straightforward signal 
of interacting pairs. 

For comparison, the GZ1 Merger class is plotted in 
Fig. [3] Since GZ2 uses a subset of the GZ1 sample, GZ1 clas- 
sifications are available for all of the galaxies in our sample 
and the same galaxies are included in Figs. [2] and [3] The GZ1 
and GZ2 Merger classifications display the same behaviour 
with decreasing projected separation, with both showing a 
dip at smallest separations (r p < 10 kpc). For GZ2 
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many of these Merger votes are exchanged for the Irregular 
and Disturbed categories, whereas for GZ1 it appears that 
the votes go to Spiral- CW and Spiral- AC W classifications, 
perhaps due to the formation of tidal arms. It is possible 
that similar effects have been present, but gone unnoticed, in 
previous work using more traditional classification schemes. 
Despite the presence of a significant signal for the control 
sample, in most separation-bins the \ow-AV pairs display 
higher Merger vote fractions than high-AV pairs, and hence 
a sensible correction for this projection bias may be applied. 

The parameters of the fitting functions, quantifying the 
amplitude and scale of the trends shown in Fig. [2] and [3] are 
provided in Table [1] 

4.2 Spiral arms and tidal tails 

Although the 'Odd' GZ2 questions were targeted directly at 
identifying interactions, there are other morphological in- 
dicators available. Tidal interactions frequently induce or 
enhance spiral arm and tidal tail features. Indeed, we find 
that a number of the GZ2 answers regarding spiral features 
display a clear dependence on pair separation. 

Fig. [4] presents the trends in 'Spiral' GZ2 morphologi- 
cal features as a function of pair separation. As in Fig. [2] 
the main panels plot the mean vote fraction, F(r p ,A), of 
AV < 500 km s _1 close pairs, after correcting for projec- 
tion bias effects using the 1000 < AV < 5000 km s _1 con- 
trol sample (see 13.2)) . This quantity represents the probabil- 
ity, P(r p ,A\M), of a galaxy in this sample displaying the 
specified morphological attribute, purely as a result of be- 
ing in a pair with separation r p , provided it is meaningful 
to ask about that feature. The raw vote fractions, f(r p ,A), 
from which the corrected quantities are determined, are also 
shown in the inset panels. Although we still see some trends 
with r p in the control sample, the projection biases are 
much less severe than those for the 'Odd' features considered 
above. 

We first consider the probability that galaxies 
will be classified as displaying spiral features, given 
that they display any features and are not edge-on 
disks, F(r p , Spiral=Yes) = P(r p , Spiral = Yes | Features D 
Edge-on=No). We find that this increases significantly with 
decreasing projected separation, beginning around r p < 100 
h?o kpc. We see that the trends in F(Spiral=Yes) and 
F{2 Arms) are similar, indicating that the increase in the 
probability of close pairs presenting spiral features is ac- 
companied by an enhancement in the probability of those 
spiral features being two- armed. This increase in the prob- 
ability of 2 Arms is at the expense of the probability that 
the number of arms cannot be discerned. Together, these 
indicate a general strengthening of two-arm spiral patterns 
in close pairs on an exponential scale of ~ 70 Ji^q kpc. This 
is consistent with observations and simulations which see an 
enhancement in spiral arm strength in interacting systems, 
often accompanied by an in crease in star formation activity 
(e.g. JSundelius et al.|[l987h . 

We can also consider the occurrence of different num- 
bers of spiral arms. At small separations (r p < 20 /if x kpc) 
the probability of a galaxy in a low AV^ pair displaying a 
single spiral arm, F(l Arm) increases sharply. There is per- 
haps a suggestion that F(r p , 2 Arms) decreases somewhat as 
1 Arm increases at small r p . The other answers to the 'How 



many spiral arms are there?' question show no significant 
change with decreasing pair separation. 

Finally we examine how the winding tightness of any 
spiral pattern depends on pair separation. In Fig. [4] we see 
that the probability of a galaxy being classified as having 
Loose Winding Arms (hereafter LWA) increases for r p < 100 
h?Q kpc. On the other hand, the Medium Winding Arms and 
Tight Winding Arms classes show no significant change with 
r p . Although it has been known for some time that strong 
tidal interactions between galaxies can produce either 1 or 
2 loosely wou nd tidal arms, depending on the orbital pa- 
rameters (e.g-, Thomasson et all 1 19891 : iHoward et al. I ll993l : 
iBarnes fe Hibb ard 2009; Barnes 201 lh . this is the first study 
to determine the observability of these features as a function 
of separation. 

The 1 Arm and LWA features appear to be particularly 
robust indicators of pair interaction, as the control sam- 
ples display opposing behaviour to the physically associated 
pairs. In the absence of a true interaction, galaxies in pro- 
jected pairs are actually less likely to be classified as having 
1 Arm or LWA. Given the limited impact of projection bias 
effects for these quantities, we could in principle use the low 
AV pair trends in arm winding and arm number directly, 
without the need for remove the r p dependence of the con- 
trol sample. However, performing this correction gives us a 
more reliable and quantitative measurement of the trends. 
A price of this correction is noise added due to the limited 
size of the control sample. 

It was originally envisaged that these questions would 
primarily provide information regarding the usual spiral 
arms, but it seems clear that in the case of close pairs they 
are revealing additional information. Observations and sim- 
ulations of strong tidal galaxy interactions frequently show 
extended, asymmetric tails (see the discussion and refer- 
ences in Section [TJ. We therefore interpret the trends in 
F(r p , 1 Arm) and F(r p , LWA) as unambiguous signatures 
of tidal tails caused by interactions between close pairs. 

We also see an increase in the Features and Edge- 
on=No answers with decreasing pair separation. The in- 
crease in Features appears to be a result of galaxies that 
would have otherwise appeared smooth and featureless gain- 
ing enhanced spiral arms or tidal features in close pairs. 
Similarly, spiral galaxies which are being tidally disturbed 
may develop warped disks, and therefore be less likely to be 
classified as Edge-on=Yes. This could have the additional 
effect of preserving the visibility of the spiral arms, despite 
the high inclination. Also, at least in some cases, it appears 
that participants may interpret a warped edge-on disc as 
LWA. 

The amplitude and scale of the trends shown in Fig. 2] 
in terms of the parameters of the best fit of Eqn. [8j are 
provided in Table [1] along with the corresponding values for 
Features and Edge-on=No. 



4.3 Barred galaxies 

The identification and properties o f barred gala x ies in 
Gala xy Z oo is studied in d etail by iMasters et al] (l201lL 
l2012h and lHovle et~aT1 faOllh . The dependence of GZ2 bars 
on environment is presented in ISkibba et all (l2Q12h . Here we 
continue our focus on the relationship between morphology 
and close pair interactions. 
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Figure 4. Mean vote fractions for the Spiral=Yes answer and answers to the subsequent 'arm number' and 'arm tightness' questions 
which show trends versus r p : 1 Arm, 2 Arms and Loose Winding Arms. The main plot in each panel shows the mean vote fraction of 
AV < 500 km s _1 close pairs, corrected for projection bias, i.e. F(r p ,A) = P(r p , A\M (points with error bars), together with best fit 
(blue line) to these points using equation [8] and the corresponding la confidence region (green lines). The insets plot f(r p , A) for the 
AV < 500 km s _1 close pairs (blue) and 1000 < AV < 5000 km s _1 control sample pairs (red), from which F(r p ,A) is calculated. 



In Fig. [5] the corrected mean vote fraction for the 
Bar=No answer, F(r p ,Bar — No), is plotted versus pro- 
jected separation. This trend was found to be rather noisy 
so we use 20 kpc r p bins to more clearly represent the 
data. There appears to be a increase in F(r p ,Bar = No) 
for r p < 20 kpc in the low-Ay pairs, while the high- 
Ay control sample pairs show no significant change with 
projected separation. 

This result indicates that bars are suppressed, rather 
than triggered, by strong tidal interactions. It is possible 
that part of the observed trend is a result of bars becom- 
ing less noticeable in interacting pairs due to increased star 
formation and looser spiral arms, but even so there can- 
not be a strong enhancement of bars in our interacting 
pairs. This is in agreement with the preliminary findings of 
iMendez-Hernandez et al.l (|2Q1 lh who compared isolated and 
paired galaxies and found the bar fraction to be ~ 43 per 
cent for isolated galaxies, but only ~ 20 per cent for pairs 
(where b ars are identified in ellipse fits to the isophotes). 
Similarly, iLee et al.1 (l2012h find the fraction of visually clas- 
sified strong bars decreases at small pair separations. 

At first this seems contrary to the expectation 
that bar features are excited by gravitational interac- 



tions, a s indi c ated by ma n y simulation studies (e.g., 
i Noguchil [l988l: IMoore et al.1 ll996l : IMoore. Lake, fc Katzl 
1998; iRomano-Dfaz et al.l I2008T L although in simulations 
including gas, tidal interactions appear less able to 
generate bars (e.g., iBerentzen et al.l l2004h . There have 
also been observational indications of bar enhancement 
in dense environments and intera cting systems (e.g., 
lElmegreen, Elmegreen, fc B ellin 1990). However, disk galax- 
ies in denser environments are more massive and possess 
redder colours and earlier morpholog ies, all of which als o 
correlate with the presence of bars (Masters et al. 1201 ih . 
In a detailed study of this issue, ISkibba et al.1 (|2012h found 
a substantial enhancement in the appearance of bars in 
denser environments, but concluded that the majority of 
the effect was attributable to the dependence of colour and 
stellar mass on envi ronment. How ever, a significant correla- 
tion remains, which ISkibba et"aD argue indicates that minor 
mergers and tidal interactions increase the appearance of 
bars by triggering disc instabilities. They also find the bar- 
environment correlation to decrease at small separations, 
becoming statistically insignificant at r p < 150 kpc. 
While not identical to our close pair results, this is nev- 
ertheless consistent. Simulations also demonstrate the sup- 
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Table 1. Best fit results to the plots of F(r p ,A) in Figs. [2] [3] 0] 
and [5] using Eqn. [8] Edge-on=No and Bar=No are fits to trends 
with 20 h^Q kpc bins, Spiral, Loose Winding Arms, 1 Arm, 2 
Arms and Other use 10 h^ kpc bins, while Odd, Merger, Irreg- 
ular, Disturbed and Features use 5 h^Q kpc bins. 
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Figure 5. Mean vote fractions for the Bar=No answer. The main 
plot shows the mean vote fraction of AV < 500 km s _1 close 
pairs, corrected for projection bias, i.e. F(r p ,A) = P(r p , A\M 
(points with error bars), together with best fit (blue line) to these 
points using equation [8] and the corresponding la confidence re- 
gion (green lines). The insets plot f(r p , A) for the AV < 500 km 
s _1 close pairs (blue) and 1000 < AV < 5000 km s _1 control 
sample pairs (red), from which F(r p ,A) is calculated. 



pression or destruction of bars in strong tidal interactions 
(Ber entzen et al.l 12003). However, most simulation work fo- 
cuses on the final result of pair interactions, rather than 
the evolution of morphological features over the course of 
the interaction. In the low-mass ratio pairs we study in this 
paper, the final result of the interaction may often be the 
destruction of any bars along with their entire disk. Theo- 
retically it seems that destruction of bars at earlier stages 
of pair interactions is feasible, and the observations, in this 
paper and by others, appear to indicate this. The emerging 
picture is therefore that moderate interactions, with high- 
mass ratios (i.e., minor mergers and harassment), promote 
bar formation, but that stronger interactions suppress the 
appearance of bars. 



Figure 6. The trend in the probability of observing Loose Wind- 
ing Arms with projected separation, F(r p , LWA), for pair sam- 
ples selected with different ranges in AV: < AV < 100 km s _1 
(black), 100 < AV < 200 km s _1 (purple), 200 < AV < 300 km 
s _1 (blue), 300 < AV < 400 km s _1 (green), and 400 < AV < 
500 km s _1 (yellow). The low AV pairs clearly show stronger 
signs of interaction than the high AV pairs. 



4.4 Dependence of tidal effects on line-of-sight 
velocity difference 

So far we have been considering the r p dependence for a 
fixed sample of pairs with AV < 500 km s _1 , versus a con- 
trol sample with 1000 < AV < 5000 km s _1 . However, 
we would also expect a relationship between physical sep- 
aration, and hence interaction strength, and AV. This is 
explored in Fig. [6] in which we plot F(r p , LWA) for pairs in 
several bins of AV. Again, we fit the data with Eqn.[8]in an 
attempt to quantify the trends. Due to the limited statis- 
tics, these plots are relatively noisy and the fit parameters 
are sometimes quite uncertain. Nevertheless, some trends 
are clearly seen in both the raw data and the fits. 

Signs of interaction appear strong for the smallest AV 
pairs, and weaken with increasing AV. We see almost no 
signs of interaction for AV > 500 km s" 1 . In Table [2] we see 
that for AV > 300 km s _1 the fitting parameter a (repre- 
senting the excess LWA likelihood at r p = /if x kpc) de- 
creases sharply. As AV decreases, there is a hint that pairs 
show signs of tidal interactions further out. This would be 
consistent with our expectations: pairs with low AV must be 
predominantly interacting in the plane of the sky, whereas 
at higher AV, pairs will be increasingly interacting along 
the line-of-sight. 

In future work we intend to treat AV and r p on similar 
terms, quantifying the dependence of observed morphologi- 
cal features as a function of both quantities. 



4.5 Identifying probable interacting galaxies 

As we have seen above, the strongest indication of in- 
teracting galaxies in physically associated pairs is an en- 
hanced probability of being classified with Loose Winding 
Arms. Previously we have studied this signal statistically, 
averaged over many galaxies in bins of r p . We now at- 
tempt to use this morphological signature to identify galax- 
ies which are likely to be interacting, by selecting those 
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AV (km s -1 ) 


a 


b (h 7 Q kpc) 


c 


< AV < 100 


0.29±0.04 


39±5 


-0.007±0.003 


100 < AV < 200 


0.18±0.04 


35±9 


-0.000±0.004 


200 < Ay < 300 


0.26±0.08 


24±7 


-0.009±0.004 


300 < Ay < 400 


0.12±0.05 


37±18 


-0.018±0.005 


400 < Ay < 500 


0.05±0.11 


21±44 


-0.009±0.004 



Table 2. Best fit results to the plots of F(r p , LWA) in Fig.[6]for 
a range of Ay bins using Eqn. [8] The fits are done using 10 
kpc bins. 



which have f(LWA) above a certain threshold. After visu- 
ally examining ~ 100 galaxies with a range of thresholds 
we found that f(LWA) > 0.6 is sufficient to reliably iden- 
tify galaxies with tidal features, provided there are at least 
two LWA votes. (A single vote may occasionally be spuri- 
ous.) Given the roughly similar number of times each object 
has been classified, the requirement of at least two LWA 
votes may be expressed as a threshold on p(M). In this 
case, p(M) = p(Features D Edge-on =No D Spiral=Yes) = 
f (Features) f(Edge-on=No)f(Spiral= Yes), which we here- 
after refer to as p(FNS). Most GZ2 objects have at least 
30 classifications, so applying a threshold of p(FNS) > 0.1 
means that the question How tight are the spiral arms? 
has received at least three answers, at least two of which 
must have indicated Loose Winding Arms if f(LWA) > 0.6. 
As described in Section 13.21 we denote this selection as 
f(LWA | p(FNS) > 0.1) > 0.6. 

The chosen thresholds are a reasonable compromise. If 
we choose lower thresholds we select potential interacting 
galaxies with more subtle tidal features, but at the same 
time increase the contamination from non-interacting pairs. 
The number of galaxies which satisfy our criteria will in- 
evitably represent only a fraction of the galaxies which are 
truly interacting. Some galaxies may possess tidal features 
that are unobservable (i.e. do not result in a significantly 
elevated f(LWA) due to their low surface brightness, an un- 
favourable sky orientation (i.e. edge-on galaxies), or other 
observational limitation. Other truly interacting galaxies 
may possess only weak, or entirely absent, morphological 
signatures of that interaction, due to t he orbital parame- 
ters o f the interaction. As discussed by iToomre & Toomrel 
galaxies which are rotating in the same sense as the 
orbital pass of the companion will form tidal arms which 
are much more pronounced than if the galaxy is rotating in 
the opposite direction. Finally, galaxies which display Loose 
Winding Arms have probably already undergone at least 
one close pass of their companion. Some close pairs will be 
in the early stages of their interaction and thus have not yet 
formed tidal arms. 

We visually examined all of the galaxies in pairs with 
r p < 200 h^Q kpc and Ay < 500 km s _1 and selected to be 
interacting with f(LWA \ p(FNS) > 0.1) > 0.6. Essentially 
all of these objects show obvious signs that they are inter- 
acting with a companion. In many cases the companion is 
also selected as interacting. However, galaxies can belong to 
multiple pairs and most of the interacting galaxies in large 
separation (r p > 100 kpc) pairs were actually found to 
have a closer companion that is likely to be the true cause 
of the interaction. In these cases the companion at larger 
separation typically shows no signs of interaction, and may 



not be physically associated. This is also the case for some 
very close pairs, where the closest companion is actually an 
interloper, and the interacting companion is at a larger pro- 
jected separation. Also, due to the redshift incompleteness 
of the SDSS, especially at small angular separations, some of 
the interacting companions of galaxies identified as having 
Loose Winding Arms will be missing from our sample. 

Table [3] shows examples of galaxies with Loose Winding 
Arms, i.e. f(LWA \ p(FNS) > 0.1) > 0.6, for a range of r p 
up to 100 h^Q kpc. We see that the high mass pairs consist 
mostly of early type galaxies, while the lower mass pairs are 
mostly mixed and late- type pairs with bluer colours. 



4.6 Accounting for interlopers in the frequency of 
galaxy pair interactions 

In this section we demonstrate how morphological indica- 
tors of interaction may be used to to refine estimates of the 
frequency of galaxy pair interactions. We employ a mass- 
limited sample, selected by imposing a minimum mass limit 
of 1O 9 5 M and redshift limits of 0.01 < z < 0.05, giving a 
sample of 44064 galaxies. The highest mass galaxies in our 
sample have M* ~ 10 1 1 5 Mq. We consider all galaxy pairs 
with Ay < 500 km s _1 and r p < 300 kpc, in 5 kpc 
bins. 

As we have seen, a significant fraction of the galaxies 
in these pairs will not be truly interacting. Previous studies 
have found this c ontamination to a s trong function of pro- 
jected separation ([Alonso et al.ll20o3 ; IPerez et al.l l2006) and 
galaxy luminosity (a nd thus mass), with th e faintest pairs 
being most affected (IPatton fc Atfieldl l2008h . We therefore 
use the morphological indicators of interactions, explored 
above, to estimate and correct for the influence of interlop- 
ers on our determination of the interaction rate. 

We calculate the average number of companions at sep- 
aration r p (and within Ay < 500 km s _1 ) for each galaxy 
in the sample, N c (r p ), following the method in Section [331 
(with u>i n t = 1)- This method accounts for the various sam- 
ple selection issues to produce a corrected estimate of the 
average number of close companions per galaxy. As some of 
these companions will not be physically interacting, this is 
an upper limit on the number of interacting companions per 
galaxy. 

We also measure the average number of companions in 
pairs displaying Loose Winding Arms, N cLW a{t v ). This is 
achieved in an identical fashion, but limited to probable in- 
teracting galaxies where at least one member of a pair is 
identified as having Loose Winding Arms, using the criteria 
developed in Section |43] i.e. w int = 1 if f(LWA | p(FNS) > 
0.1) > 0.6 for either member of the pair and Wi n t = other- 
wise. 

The fraction of companion galaxies which show indi- 
cations of possible interaction (i.e. with f(LWA \p(FNS) > 
0.1) > 0.6), can be calculated for each r p bin using, 



J^LWA^p) 



N c lwa(t v ) 
Nc(r P ) 



(12) 



This fraction is related to the true interacting fraction, 
^intt^p), such that 



^int(j p ) 



FLWA{r p ) — Tlwa(t p — > oo) 

-Pintobs 



(13) 
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Table 3. Galaxies selected to be interacting members of a pair, using f(LWA \p(FNS) > 0.1) > 0.6 and AV < 500 km s -1 . The figures 
along the top indicate the approximate r p (in kpc) of the pairs in each column. The other figures give the stellar masses (in Mq) 
of galaxies in the adjacent image. The top and bottom rows present examples with high and low stellar mass, respectively. The width of 
these figures is ~ 110 kpc, based on the average redshift of the pair. 



where Tlwa(t p — > oo) gives the fraction of galaxies passing 
our LWA selection in the absence of pair interactions probed 
by our sample selection. These galaxies are mostly interlop- 
ers, associated in projection with truly interacting pairs: a 
galaxy with LWA interacting with one physically close com- 
panion could also have additional, non-interacting, compan- 
ions at any r p . These companions would be (falsely) counted 
as interacting in Eqn. [12] Some galaxies may also display 
LWA features due to intrinsic properties of the galaxies or 
due to interactions with companions possessing masses lower 
than our sample selection. Both these cases are accounted 
for by subtracting Tlwa(t p — > oo) in Eqn. 1131 

The factor Pintobs converts the fraction of galaxies in 
pairs displaying LWA into the fraction of truly interact- 
ing galaxies. It is the average probability of a true physi- 
cal interaction resulting in an observable LWA signature in 
our dataset. In principle, Pmtobs could be a function of r p , 
although we expect it to vary slowly. Remember that our 
working definition of an 'interaction', from Section [1] is that 
a galaxy has experienced a significant tidal force, compared 
to its gravitational binding force, averaged over the previous 
dynamical time. For a given dataset, the level of tidal force 
that is deemed 'significant' is that which results in observ- 
able features in the most favourable circumstances. Varia- 
tions of Pintobs from a constant with r p are therefore only 
expected from secondary effects. However, this assumption 
would greatly benefit from being tested with simulations, 
which would potentially result in a refined functional form 
for Pintobs- For the time being, we assume a constant Pintobs- 
Although the LWA class is a rather indirect indicator of tidal 
tails, Pintobs accounts for the difference between the number 
of objects actually counted and the number missed. Future 
refinements to Pintobs could included a more careful consid- 
eration of the conditions under which interactions produce 
tidal tails that would be classified as loose spirals. 

Figure [7] plots Tlwa{t p ) and a fit to the trend using 
equation [8] The best fit gives parameter values olt lwa — 
0.236, bjr LWA = 20.145 kpc and cj? LWA = 0.035. At large 
r p the curve levels off to a constant value, Tlwa(t p —> oo) = 



c ^lwa - Therefore, 

Jint(rp) = a ^wA^M-r P /bT LWA ) (M) 

-Pintobs 

We may place strong constraints on Pintobs by noticing 
that .Fmt(r p ) must lie on the interval [0, 1]. Therefore, a con- 
stant Pintobs must lie on the interval [clt lwa , 1] , and hence 

^int(rp) = [ajr LWA ,l] • exp(-r p /bjr LWA ) , (15) 

where the square brackets denote an interval. This provides 
the means to calculate strong limits on the true frequency 
of interactions from our observations of galaxies with LWA 
features. 

We can determine an estimate of the average number of 
interacting companions for galaxies in our sample, Af cint (r p ), 
by applying an r p -dependent weight, it; m t, in the method- 
ology of Section 13.31 This weight is simply the fraction of 
companions at r p that are truly interacting galaxies, so 
Wint = ^int(r p ). For our sample, using the parameters of 
the fit to Tlwa^p) in Fig. we derive the interval, 

Mint = [0.236, 1.0] • exp(-r p /20.145 kpc) . (16) 

Usually when one calculates the close companion fre- 
quency, one must apply an r p limit. Using u>i n t, however, 
weights companions by the likelihood that pairs with that 
projected separation are truly interacting. One can therefore 
integrate to large r p and determine a total frequency of in- 
teracting companions. To illustrate, without the interaction 
weighting (i.e. Wi n t = 1 in Eqn. [9]), considering pairs with 
small projected separations gives N c (r p < 30 h^Q kpc) = 
0.028 ± 0.002. Including the w int given by Eqn. [T6l results 
in N cint (r p < 30 kpc) = [0.0033 ± 0.0010, 0.014 ± 0.004], 
where the interval brackets the allowed range of probabil- 
ities, Pintobs, that an interaction results in an observable 
LWA feature. Comparing Af c and the upper limit for Af cint 
implies that at least 49 =b 14 per cent of pairs with r p < 30 
/ifo 1 kpc, AV < 500 km s -1 and M* > 1O 9 5 M , are 
non-interacting int erlopers. This agrees with the results of 
IPatton et aD d2000h , which found using visual classifications 
that ~ 50 per cent of galaxies were interlopers for a similar 
range of r p , AV and M*. 
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Figure 7. Pairs containing at least one member with 
f(LWA\p(FNS) > 0.1) > 0.6 divided by the total number of 
pairs in each r p bin. Results are given for AV < 500 km s _1 
(black line), together with a fit using equation [8] (blue line) and 
the corresponding la confidence region (red lines). 

If we select all pairs up to our limiting projected sep- 
aration, we find that N c (r p < 300 kpc) = 0.62 ± 0.004, 
i.e. most galaxies have a companion within this distance. 
On the other hand, the average number of physically inter- 
acting companions per galaxy is N C i nt (r p < 300 Ji?q kpc) = 
[0.0048 ± 0.0014, 0.021 ± 0.006]. This implies that > 96 per 
cent of pairs with r p < 300 h^Q kpc, AV < 500 km s _1 and 
M* > 1O 9 5 M are not actually interacting. Our method 
enables one to estimate the fraction of interacting galaxies 
without requiring an arbitrary cut-off in projected separa- 
tion, and without any need for further contamination cor- 
rections. 

Note that while we have developed a more sophisticated 
treatment of r p trends here, we are still using a simple cut in 
AV. In principle, the method outlined here could possibly be 
extended to determine an interaction weight, w-mt, with de- 
pendence on both r p and AV^. There is also the potential of 
further constraining the interaction frequency by combining 
multiple indicators of interaction, rather than LWA alone. 
This approach could perhaps also enable the identification 
of pairs at different stages in their interaction, or with spe- 
cific orbital characteristics. We leave all of these challenges 
for future work. A paper using the methods outlined above, 
to explore the dependence of the interaction frequency on 
stellar mass and environment, is in preparation. 



5 SUMMARY AND DISCUSSION 

In this paper we have examined a variety of morphological 
signatures of interaction between galaxy pairs, and demon- 
strated how the trends in these observable features, as a 
function of projected separation, can provide a refined es- 
timate of the frequency of pair interactions in the galaxy 
population. We consider an 'interacting' galaxy to be one 
which has experienced a significant tidal force, compared 
to its gravitational binding force, averaged over the previ- 
ous dynamical time (see Section [1]). The tidal force deemed 
'significant' depends upon the properties of a given observa- 
tional dataset. 



We began by presenting our sample and the methods 
we employ, and particularly, in Section 13.21 discussing the 
information provided by Galaxy Zoo 2 and its interpretation 
in terms of the probability that a given galaxy is observed to 
possess a particular set of morphological features. We also 
presented a method to correct for 'projection bias', an effect 
whereby the signal of certain morphological features may 
depend on the apparent separation of galaxy pairs, even in 
the absence of any possible physical associations between 
pair members. 

In Section 14.11 we considered questions from GZ2 de- 
signed to identify Odd features, including answers, such as 
Merger, which were intended to identify signs of interaction. 
We found that these classifications suffer from a strong pro- 
jection bias. For example, galaxies with small projected sep- 
aration, but very large velocity offsets, tend to have a spu- 
riously large Merger signal. Furthermore, the way in which 
the question for these Odd features was arranged, allowing 
only one of the available options to be selected, results in 
cross-talk between the Odd categories, in terms of both true 
signal and projection bias, which complicates their inter- 
pretation. We find that this projection bias is also present 
for the GZ1 Merger classification, and shows a behaviour 
very similar to its GZ2 equivalent. Previous studies which 
used the GZ1 Me rger class to identify merger candidates 
([Darg et al. 2010a : b) will have suffered to some degree from 
this issue, but the effect is probably relatively small due to 
their use of vote fraction thresholds and the fact that the 
low- Ay galaxy pairs have a larger mean Merger vote frac- 
tion at most projected separations relative to the control 
sample pairs. In future iterations of Galaxy Zoo, and other 
visual classification efforts, it would be preferable to keep 
questions regarding different types of features distinct, or 
allow multiple answers to be selected for a single question 
when the relevant features are not mutually exclusive. 

Nevertheless, the Odd classifications do provide useful 
information on the reality of galaxy interactions, particu- 
larly once a correction for the projection bias is applied by 
reference to a control sample of pa i rs with large velocity off- 
sets. As discussed by iDarg et al.1 (j2010aT ) with relation to 
GZ1, the GZ2 Merger class (and also Odd=Yes) primarily 
selects interacting galaxies at small projected separations. 
It therefore mainly probes close passes and the later stages 
of mergers. We further find that the Irregular and Disturbed 
classes can identify interacting galaxies with very small pro- 
jected separations, which may be either at an advanced stage 
of merging or aligned along the line- of- sight. 

We have searched all the GZ2 classifications for trends 
with projected separation, and find significant signals with 
respect to the presence and form of spiral arms. The ob- 
servability of spiral arms (Spiral=Yes), and particularly the 
dominant 2 Arms cl ass, is enhanced for close pairs on a scale 
of ^ 70 kpc. [Darg et al.1 feOlOsJ ) find that the spiral- 
to-elliptical ratio for galaxies classified as mer gers in GZ1 
is appr oximately twice the global ratio, and in IDarg et all 
(| 201 Obi ) conclude that this is due to the longer time-scale 
over which spiral mergers are detectable compared to ellip- 
tical mergers. Our results show that this ratio can also be at 
least partly explained by the enhancement and formation of 
spiral arms in interacting galaxies. More unusual spiral arm 
features also present a trend with r p . The occurrence of One 
Arm spirals dramatically increases for small separations (on 
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a scale of r p < 20 kpc), while Loose Winding Arms 
show the strongest increase (operating on an intermediate 
scale of r p < 30 kpc). 

There are two principal ways in which spiral- like fea- 
tures can be created through tidal interactions. Tidal per- 
turbations can instigate or amplify instabilities in gas disks, 
leading to the formation or enhancement of star formation 
in spiral ar ms similar to th ose seen in isolated galaxies (e.g., 
as seen bv lXu et alJl2010h . Tidal forces can also strip stars 
and gas out of the galaxies, forming tidal tails, counter tails 
and bridges , which may or may not harbour star- format ion 
(e.g., iMullan et al.ll201lh . We appear to detect both signa- 
tures: an enhancement of 'normal' spiral arm features, oc- 
curring at large projected separations, with the signatures of 
stronger tidal interactions becoming increasingly prevalent 
at smaller separations. The tidal nature of the Loose Wind- 
ing Arms features is confirmed by examination of typical 
images, such as those in Table [3] It is clear that many of the 
galaxies which have significant probability of Loose Wind- 
ing Arms, especially those with higher stellar mass, are red, 
early- type galaxies. This indicates that the loose spiral fea- 
tures that are observed in these galaxies, and probably also 
those same features in star- forming, late- type galaxies, are 
the result of tidal stripping. In galaxies with sufficient cold 
gas, there will almost certainly be star formation in these 
tidal spiral features, and indeed in Table [3] we see that sev- 
eral of the gal axies possess very blue loose spiral arms. Sim- 
ulations (e.g ., iToomre fe Toomrelll972l : iHoward et all 1 1993: 
Barnes 2011) indicate that tidal features, such as those ap- 
parently identified by the GZ2 Loose Winding Arms class, 
are indicative of a stage between close passes, primarily be- 
tween 1st and 2nd pass, when pairs can still attain rela- 
tively large separations. The Loose Winding Arms features 
are thus probing the early stage of mergers and pair inter- 
actions. 

The onset and appearance of tidal arms are known 
to depend on the geometry of the encounter. Numeri- 
cal studies demonstrate that in-plane, prograde encoun- 
ters produce the most symmetrical two-sided disturbances, 
while polar encounters give the most one-sided distur- 
bances, and re trograde encounters are the last to make 
tidal tails (e.g.. [Thomasson et al.lll989l : IHoward et al.l fl993: 
iBarnes fc Hibbardll2009l ; lBarnedl201lh . Retrograde encoun- 
ters also produce the gre atest increase in star formation 
efficiency (|Cox et a l. 2008). Considering this, the galaxies 
which are selected as having 2 arms and Loose Winding 
Arms in GZ2 are likely the result of prograde, in-plane en- 
counters, while galaxies identified as displaying 1 arm are 
likely the result of polar or retrograde encounters. Our ob- 
served separation scales for these different features are con- 
sistent with this interpretation (see Table [J). 

When comparing our results to other studies which look 
at the onset of tidally induced changes in interacting galax- 
ies, we find that answers to Galaxy Zoo questions regarding 
spiral arms detect changes at separations similar to stud- 
ies of tidally induced star formation, as discussed in Section 
[TJ The Loose Winding Arm class begins to detect interact- 
ing galaxies around r p < 120 /ifg 1 kpc, which is similar to 
the separation scale associated with induced star form ation 
(iNikolic et al.ll2004l : uTet al.ll2008l : IPatton et alll201lh . The 
star formation detected at these large separations is rela- 
tively weak, while a strong increase is observed for r p < 40 



k pc (iLi et al.1 2008 |l Ellison et al.l 120081 : iRobaina et al.l 
120091 : IPatton et aD " 



120111 ) 



This corresponds to scale for 
which we observe an enhancement of the Merger class. 
Quantitative morphological meas urements also typically 
present signals on these scales (iHernandez- Toledo et all 
120051 : iDe Propris et al.ll2007l : lEUison et alJl20ld ). Logically, 
kinematic disturbance must precede morphological distur- 
bance and so it might seem that star formation, if it 
is triggered by kinematical perturbations, should be an 
earlier indication of interaction than morphology (e.g., 
iBvrd fe Howa rd 1992). However, this paper shows that some 
morphological signatures are as sensitive as enhanced star- 
formation, and more unambiguously related to interaction. 

In Section T4. 31 we found that the likelihood of a bar be- 
ing observed decreases sharply for pairs with projected sep- 
arations r p < 20 h^Q kpc. Bars are t hought to be created 
throu gh periodic orbital resonance ([Bournaud fe Combed 
l20Q2h and are known to initiate radial gas inflows, which 
in th e end act to destroy or weaken the bar struc- 



see Sellwood & Sanchez 


(120101). or 


ISellwood & Wilkinson! ( 


1993h). Gas 



inflows, perhaps to- 
gether with the enhancement of bar features, caused by 
tidal perturbations in the early stages of major interac- 
tions may similarly act to rapidly destroy any pre-existing or 
trans ient bars (e.g., iBerentzen et al. 2003; D i Matteo et al.l 
120071 ). Our results suggest that this is indeed the case, 
with the appearance of bars being strongly suppressed 
in close pairs, in agreement wit h oth er recent studie s 
by iMendez-Hernandez et al.l (|201 lh and iLee et al ] (|2012h . 
Eventually, the violent reorganization of stellar orbits in the 
later stages of many major interactions (i.e. mergers) must 
act to erase any orbital resonances which created the bar. 

In Section [4. 5 1 we focus on using the presence of Loose 
Winding Arms to identify probable interacting galaxies. 
These criteria are then used in Section 14.61 to constrain 
the frequency of galaxy pair interactions, without requir- 
ing an arbitrary cut-off in projected separation or any fur- 
ther corrections for contamination of our close pair sample. 
We find that the fraction of galaxies with M* > 1O 9 5 M 
and 0.01 < z < 0.05 that are in truly interacting pairs with 
AV < 500 km s -1 is in the range 0.5 ±0.1-2.1 ±0.6 per cent. 
The limits correspond to assuming the maximum and mini- 
mum permitted probability, Pmtobs, that interacting galaxies 
produce observable Loose Winding Arms features, respec- 
tively. We expect simple extensions of our technique to lead 
to significantly tighter confidence intervals in future work. 

It is difficult to precisely compare our estimate of the in- 
teracting galaxy fraction to other studies, due to the range of 
different methods employed. Althou gh the close pair fractio n 
is mostly constant with luminosity (|Patton fe At fieldll 20081 ). 
the limiting mass-ratio and projected separation, varying 
definitions for selecting pairs, and many other subtleties, 
make comparisons complicated. Given that the estimate in 
this paper is derived from a relatively simple demonstra- 
tion of combining close pair and morphological informa- 
tion, we defer such involved comparisons to future work. 
Nevertheless, we note that the major interaction fractions 
quoted by other recen t studies: e. g., 1 .1 ± 0.5, 2.1 ± 0.1, 
1.6 ± 0. 1, 1.3 ± 0.1 by iBell et al.l (l2006h: IPatton fe Atfieldl 
(120081 ): iDomingue et al.1 d2Q09h : Ixu et al.1 (l2012h . respect 
tively, are neatly bracketed by our estimate of 0.4 - 2.7 per 
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cent. This lends support for our use of GZ2 Loose Wind- 
ing Arms as indicators of pair interactions, and encourages 
confidence in the method described in Section 14.61 and the 
various assumptions we have made. 
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